Spin-orbit corrections of order ma 6 to the fine structure of (37, 35) state in 4 He + p atom. 
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Precise numerical calculation of radiofrequency intervals between hyperfine sublevels of the (37, 35) 
state of the antiprotonic helium-4 atom is presented. Theoretical consideration includes the QED 
corrections of order ma 6 to the electron spin-orbit interaction. The effective Hamiltonian is derived 
using the formalism of the nonrelativistic quantum electrodynamics (NRQED). 

PACS numbers: 36.10.-k, 31.15.A-.31.30. J- 



I. INTRODUCTION 

The high precision spectroscopy measurement of the 
hyperfine structure of the antiprotonic helium has two- 
fold interest. First, it is expected that it may be a way 
to obtain improved value of the magnetic moment of an 
antiproton. The other point is that it can be a good 
benchmark for testing QED theoretical methods for the 
Coulomb three-body bound states to a high precision. 

At present several theoretical calculations for the hy- 
perfine structure of the (37, 35) state of the 4 He + p atom 
have been performed [H, 0, HI 0] ■ Since all the results were 
obtained within the frames of the same Breit-Pauli ap- 
proximation, the major difference in the obtained data 
was either due to numerical inaccuracy of the nonrela- 
tivistic solution or due to a difference in the choice of the 
physical constants. Still they were in a good agreement 
(within the error bars of the theoretical approximation) 
with the first experimental observation of the 13 GHz 
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FIG. 1: Schematic diagram of hyperfine sublevels of the 
(37, 35) state of 4 He + p atom. 



intervals 

Recently, the ASACUSA experiment has obtained new 
precise values for the two RF transitions of the (37, 35) 
state in the 4 He + p atom [|| 0] (the notation is shown on 
the schematic diagram of the (37, 35) state on Fig. HJ): 



t+ = 12 896.641(63) MHz, 
r_ = 12 924.461(63) MHz. 



(1) 



The difference between r+ and r_ is nearly proportional 
to the antiprotonic magnetic moment and has a value 



At = 27.825(33) MHz. 



(2) 



That should be compared with the theoretical predic- 
tion 0: 



t+ = 12 896.35(69) MHz, 
r_ = 12 924.24(69) MHz. 



(3) 
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It is seen that the experimental error is more than an 
order of magnitude smaller and transitions have some 
systematic shift toward larger values. The Breit-Pauli 
approximation used so far is limited by the uncertainty 
of order 0(a 2 ) and higher order corrections should be in- 
cluded into consideration to achieve the similar accuracy 
as in the experiment. 

The effective Hamiltonian of the hyperfine interaction 
may be written (see details in [3j) 

H cS = E 1 (s e -L) + E 2 (s p - • L) + E 3 (s e • s p ) 

+E i {2L(L + l)(s e -s p ) (4) 

-3[(Sp-L)(s e -L) + (s e -L)(s p -L)] }. 

To get improved values for the r + and r_ transitions 
one needs to get contributions of the next to the leading 
order for the electron spin-orbit interaction coefficient E\ . 
It may be done within the framework of the NRQED 
formalism . Some details of the derivation of the ma 6 
order contributions may be found in 0]. 
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II. CORRECTIONS OF ORDER ma b TO THE 
EFFECTIVE HAMILTONIAN OF THE FINE 
STRUCTURE. 

In what follows we use the notation: Ri, R 2 , r e , and 
Pi, P2, p e are coordinates and impulses of three particles 
in the center of mass frame, where subscript 1 stands for 
a helium nucleus and 2 for an antiproton. We also make 
use Ti = r e R; , i = 1,2 for coordinates of an electron 
with respect to one of the nuclei. 

According to the NRQED the effective Hamiltonian 
includes the three new interactions, which contribute to 
the electron spin-orbit term at the ma 6 order. Two are 
relativistic corrections of order v 2 /c 2 to the vertex func- 
tions for the spin-orbit and Fermi interactions (see Q) 
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(5) 



The third is the seagull vertex interaction with one 
Coulomb and one transverse photon lines: 



V, = e 4 




(6) 



+ (1 «-» 2) 



the sources may be two different nuclei or may coincide. 

To get a complete set of corrections one needs to take 
into account the second order contribution as well 
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where 



2m 2 



[Z 1 6(r 1 ) + Z 2 S(r 2 ) 



(7) 



(8) 



Radiative corrections (form factors of the electron) 
have been already included into consideration as contri- 
butions to the anomalous magnetic moment. 

Transforming potentials Vi to the coordinate space and 
atomic units one gets: 



and 
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(9) 
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III. VARIATIONAL WAVE FUNCTION 

For numerical calculations we use the exponential vari- 
ational expansion, which has been discussed in details in 
ficj |. Namely, the wave function for a state with a total 
orbital angular momentum L and of a total spatial parity 
7r = (— 1) L is expanded as follows: 

*(R,ri)= Yl yLMi^h)Gfj 2 (R, ri ,r 2 ), 

h+l 2 =L 

N 



Gt^{R,r x ,r 2 ) 



n=l 



+D n Im[e 



where R is a position vector of an antiproton and ri is a 
position vector of an electron with respect to a nucleus; 
parameters a n are complex, and /3 n , 7„ are real, they are 
generated in a pseudorandom way, 



Re(a„) = [[\n{n + l)y/p^\ (A 2 - A x ) + A t ] 



Im(a n ) 



(12) 



in(n+l)v^| (A> 2 -A< 1 )+ A\ 

n = [[\n{n + (B 2 - B ± ) + B x ] , 

7 „ = [\_\n{n + ljy^j (C 2 - C x ) + C x ] . 

Here [^J denotes fractional part of x, and p a , pp or p 1 
are some prime numbers and Ai, Bi, Ci are variational 
parameters. 

For the initial wave function of the bound state we 
use the triple basis set with the total number of terms 
N = 2200 in expansion (jTTJ) and full optimization of 
variational parameters. That yields the non-relativistic 
energy for this state 

£ nr (37, 35) = -2.899 282 183 295 31(1) au 

The CODATA06 recommended values [Tl| have 
been adopted for calculations: nip — m p = 
1836.152 672 47(80)m e , m a = 7294.299 5365(31)m e 
and R^c = 3.289 841 960 361(22) x 10 6 MHz. 

For the intermediate states of the second order itera- 
tion the similar variational expansion (jTTJ) with various 
basis lengths TV = 520-^960 has been used. 



IV. REDUCE A SINGULARITY IN THE 
SECOND ORDER CONTRIBUTION 

The H B operator in the second order term ([7]) is too 
singular. It requires careful consideration because inter- 
mediate states should include functions with asymptotic 
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set 


[Ai,A 2 ] 


[A[,A' 2 ] 


[Bi,B 2 ] 


[Ci,C 2 ] 


1-st set 


[66.6, 87.6] 


[0.4, 5.2] 


[0.00, 2.05] 


[0.00, 0.87] 


2-nd set 


[66.0, 75.4] 


[0.0, 5.4] 


[0.94, 5.70] 


[0.00, 1.94] 


3-d set 


[66.0, 75.4] 


[0.0, 5.4] 


[5.00, 80.0] 


[0.00, 0.10] 


4-th set 


[66.0, 75.4] 


[0.0, 5.4] 


[0.00, 0.20] 


[2.00, 70.0] 


5-th set 


[66.0, 75.4] 


[0.0, 5.4] 


[90., 1000.] 


[0.00, 0.10] 


6-th set 


[66.0, 75.4] 


[0.0, 5.4] 


[0.00, 0.10] 


[80.0, 800.] 


7-th set 


[66.0, 75.4] 


[0.0, 5.4] 


[10 3 , 10 4 ] 


[0.00, 0.10] 


8-th set 


[66.0, 75.4] 


[0.0, 5.4] 


[0.00, 0.10] 


[800., 10 4 ] 



Thus the second order term may be rewritten as follows 



(H B \Q(E -H )-iQ\H so ) 



H' B \Q{E -H )- l Q\H, 



so 



(16) 



+ (UH so )-(U) (H so )). 
Matrix elements of H' B may be obtained directly from 
Eq. (|13p . Additional term to the effective Hamiltonian is 
expressed 



TABLE I: Variational parameters for eight basis sets used in 
the second order contribution calculations. 



behaviour at small distances like ~l/ri (or I/V2). The 
usual regular trial functions would result in a very slow 
convergence of AEa ■ 

In order to smooth the perturbation and to reduce the 
singularity of the intermediate wave function we may use 
transformation 



H' B =H B - (E - H )U - U(E - H ) 



(13) 



The delta-function singularity in \H b ^q) has the follow- 
ing structure 



1 

+Z 2 





-I) 







7r£(ri 



7r5(r 2 ) 



(14) 



where 1/fXi = l/m e + l/Mi. 

It is natural to take U in the form U = c\jr\ + Cijri. 
The coefficients c, may be obtained by substituting U 
into the initial Schrodinger equation 



(E Q -H )[^ + ^ 
ri r 2 



2ci 
Mi 



=£1 

M2 



7r<5(r 2 ) 



then comparing the latter expression with Eq. (fl~4")) one 

gets: 



Hi{2ni-m e ) „ 
c i = ; — 5 A, 



C2 = 



Ami 
^ 2 ( 2^2-w e ) 
4mJ 



(15) 



H^) = {UH so )-{U) (H so ) 



-(2L + S*)(H so/ 

n r 2 



(17) 



For the numerical evaluation of the second order term 
from Eq. p6[) we use the eight basis sets, where the first 
two approximate the regular part of the intermediate so- 
lution, and the remaining six sets with growing exponents 
are introduced to reproduce behaviour of the type ln(ri) 
(or ln(r 2 )) at small values of r\ or r 2 . The particular 
variational parameters used are presented in Table |TJ 



NUMERICAL CALCULATIONS 



Matrix elements in Eqs. (|g|)- (fTUa|) and ([T5 ]) -([T7 |) for 

the basis functions (jTTJ) of the exponential variational 
expansion were evaluated analytically using the recur- 
rences derived in [l§] with some modifications, which 
allowed to improve stability. The generating functions 
r_ 4! o,o(a, /3, 7), r_ 3i _i !0 (a, f3, 7) were taken from 
What corresponds to the cut-off regularization of the in- 
tegrals at rj.,2 = p, where p <C 1. 

For all vector operators the reduced matrix element is 
assumed. Here we present numerical values of some the 
most complicate operators. 



2 i r l X Pe 



Pe 



[n x Pi] 

r 3 



1.5575929 



= 3272.7020 



pi ^ * Pe ^ ) : 2.24. r )!)243 



p 2 e ^ 2 ^3 P ^ ) = -30S0.3N7!) 



(18) 
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III 


11 2 


f ^3 


[ri x p e ] 


[ri x Pi] 


[r 2 x p e ] 


[r 2 x P a ] 


r? 


rf 


r| 




20 


20 





0.2883781 


680.0299 


0.5135404 


-1204.595 


20 


20 


20 


0.2411901 


624.7925 


0.4875501 


-1177.070 


40 


20 


20 


0.2242844 


573.8206 


0.4899458 


-1192.019 


60 


20 


20 


0.2536075 


633.2949 


0.5022857 


-1194.613 


60 


40 


20 


0.2232982 


526.5736 


0.4881029 


-1187.319 


60 


40 


40 


0.2847268 


714.1212 


0.4849738 


-1160.058 


80 


40 


40 


0.2812069 


716.9230 


0.4696484 


-1159.465 



TABLE II: Convergence of the second order contribution matrix elements for the spin-orbit interaction. 



1 [ri x P 2 ] 



205.83272 



r 2 

3 I , (riPi) ) = 551.65420 



1 [r 2 x Pi] 



1391.5321 



[ri 3 X 3 r2] (r 2 P,i ) :>:, 1.4:«2.s 



(19) 



H' 



H' 



Q(E -H )- 1 Q 



[ri x p e ] 



[ri x Pi] 



0.2812 



717. 



H' 



H' 



Q(E -H y 1 Q 



[r 2 x p e ] 



[r 2 x P 2 ] 



0.4696 



-1160. 



(20) 



It is worthy to say that the second order iteration, 
even after reduction of the singularity, still reveals slow 
convergence. In Table [TT] we present results of numerical 
calculations for various sets of basis functions. The re- 
sults depend very little on increase of the 1st and 2nd ba- 
sis sets (see Table [XJ) , which represent regular behaviour. 
The following notation has been used in Table [Hi n\ is 
a number of basis functions for 3d and 4th sets, n 2 is 
for 5th and 6th sets, etc. As is seen from the Table, no 
more than two digits may be accepted with confidence as 
convergent. However, the increase of the basis sets leads 
to numerical instability, which we attribute to very large 
angular momentum of the state [L = 35) what makes 
the recursion used for analytic evaluation of the matrix 
elements to be too long and unstable for large exponents. 
The octuple precision has been used in these calculations 
and still it was not enough to provide necessary stability. 



VI. RESULTS AND DISCUSSION 

Summing up the contributions from Eq. J7]), (JSJ), and 
(I10ap - (|10b[) we obtain the ma 6 order contribution to the 
electron spin-orbit interaction: 

A.E1 = -0.000030(4) • 10~ 7 au (21) 



Thus a new value for the E\ coefficient would be 

E x = -0.552 563(4) • 10~ 7 au (22) 

where the uncertainty is primarily due to slow conver- 
gence of the second order iteration. 

Using this new value for the E\ coefficient and keeping 
E2-E4 as in [3J] one may solve the effective Hamiltonian 
dU) and get updated theoretical values for transition fre- 
quencies: 

t+ = 12897.0(1)(3) MHz 

r_ = 12924.9(1) (3) MHz (23) 
At = 27.897(0)(3) MHz 

The first error indicates the numerical uncertainty of 
present calculations, while the second one is an estimate 
of the theoretical uncertainty due to yet uncalculated 
higher order terms. As it should be At does not change 
its value comparing with previous calculations [3]. The 
values of r + and r_ leap over the experimental result and 
become to overestimate ones if we take the numerical er- 
ror as a measure of uncertainty. Still in theory we need 
to include into consideration effects of the next order in 
a, which contain terms of order (a 3 \na)Ei and are of 
the magnitude of the discrepancy. 

It is worthy to note here that the obtained value of 
AE\ is unexpectedly small. That explains rather good 
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agreement of the experiment with the results of the Breit- 
Pauli approximation. 

In order to get the improved value for At one needs 
to perform a complete calculation of all the contribu- 
tions of order ma e (m/M), which provides corrections to 
the remaining coefficients, £2 — E±, in the effective HFS 
Hamiltonian ([J). This work is in progress now. 
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